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EFFECT OF OSCILLATOR INSTABILITY ON TELEMETRY SIGNALS 


Terry L. Grant and Randolph L. Cramer 
Ames Research Center 


SUMMARY 


This paper discusses short-term oscillator instability and its effect on 
telemetry signals. Instability measurements are related to the total phase 
jitter expected in a telemetry receiver and the resulting degradation of telem- 
etry performance. The analysis of phase jitter present in the telemetry 
system of Pioneer VI through IX serves as an example. The instability is 
described by the measurement of oscillator phase jitter, and the function and 
operation of the measuring instrument are explained. 


INTRODUCTION 


In spacecraft telemetry the trend has been to use a pulse code 
modulation - phase shift keyed - phase modulated (PCM-PSK-PM) system because, 
with coherent carrier detection, such a system can theoretically provide the 
lowest error rate for a given signal-to-noise ratio of all the known modula- 
tion techniques. The key factor in obtaining this performance, however, is 
the coherent or phase-tracking detection process, and deviation from this 
ideal model quickly causes a loss of performance. It was the quest to make 
the Pioneer telemetry system perform as well as possible which motivated the 
measurement of oscillator instability and the analysis of its effect on a 
telemetry system. (Performance here is understood as the channel error rate 
versus a normalized signal-to-noise ratio.) Lack of coherence in carrier 
tracking also affects the resolution of any range and range rate measurements 
that may be made on the telemetry carrier. 

In the design of an optimum coherent, or phase-tracking telemetry system, 
several sources of phase jitter must be considered. The most obvious source 
is the additive noise in the bandwidth of the tracking loop. A second source 
of phase jitter is oscillator instability in the transmitter reference oscil- 
lator or the receiver voltage controlled oscillator (VCO) . Other sources are 
propagation anomalies that affect the phase of the carrier. The latter are 
caused by small time varying changes in the dielectric medium (ref. 1). Only 
oscillator instability as a source of phase jitter will be treated here and 
its effects on the performance of a one-way (spacecraft-to-ground) telemetry 
system will be analyzed. It then will be shown how this source of phase 
jitter forms an important bound on telemetry performance. 



This paper is a tutorial description of oscillator instability and its 
effects on a phase lock loop (PLL) . Measured phase jitter is used to analyze 
statistically the conditions under which this jitter will degrade the perfor- 
mance of a one-way (spacecraft-to-ground) telemetry system. The appendixes 
discuss the concepts of probability and statistics which are key to the analy- 
sis. Finally, the oscillator instability measuring instrument is described 
and the results of measurements on the oscillator circuits of the Pioneer 
transmitter are presented. These data formed the basis for the assumptions of 
the character of oscillator instability. (The Pioneer telemetry system 
provided the example that motivated the analysis of the effect of oscillator 
instability on a telemetry system.) 


ANALYSIS 


The key factor in obtaining optimum performance with a PSK communication 
system is maintaining a coherent reference to compare with the transmitted 
signal. Although theoretically no signal power is allowed for this reference, 
it is necessary in practice. Theoretically, both the transmitter and the 
receiver have perfect time information available; therefore nothing would be 
gained by including a time reference with the signal. In reality, sufficiently 
accurate time information is not available to allow coherent phase detection. 
The following example illustrates this point. Suppose the clocks at a trans- 
mitter and receiver measure time to a relative accuracy of ±1 ys (a difficult 
feat within the present state of the art) . Suppose further that the known 
signal is PSK modulated of the form, sin[m c (t + t 0 ) ±tt/2 ] , where t 0 corresponds 
to a known time of a maximum modulated signal. The receiver would multiply the 
incoming signal by cos w c (t + t 0 ±At) and average the result over one symbol 
period. The decision that the transmitted tr/2 rad phase shift is plus or minus 
depends on whether the average is positive or negative. It is obvious from 
figure 1(b) that the average result is zero and no information is received 




The product of I and 2 



(a) Correct time information. 


( I ) sin [* c (t + t 0 ) + tt/2] 



(2) cos w c (t + t 0 + At), At = tt/2oj c 



The product of I and 2 
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(b) Time error At = 7 t/2co c . 


Figure 1 Effect of incorrect time on phase detection. 
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if At = tt/2o) c . It can be easily shown that if At is slightly greater than 
tt/2u) c , the wrong decision would be made. Since At = 10 6 implies that 
w c < tt/2x10 6 , then a) c /27r = f c would have to be less than 1/4 MHz for coher- 
ent detection, a frequency much too low for spacecraft communication. There- 
fore the relative time must be derived from the signal itself, which is the 
function of a phase tracking receiver. In the design of such a receiver 
oscillators are not as ideal as postulated for the above example. An oscil- 
lator has phase instability, and this instability can be modeled in terms of 
the above discussion by letting At be a random process instead of a fixed 
uncertainty. Thus, even with the phase tracking receiver, the question to be 
answered is: For a given signal frequency, what is the maximum variance of 
At that allows the detection process to be considered coherent? More specif- 
ically: How does performance differ from that of coherent detection as a func- 
tion of At vs. l/f c or, equivalently, the ratio At/T c (where T c is the 
period corresponding to f c ) ? This ratio relates jitter in time to jitter in 
fractions of a cycle, but more conventionally , the ratio is multiplied by 
2tt to relate time jitter to radians of "phase" jitter. Note that frequency 
is a rate of change of phase; that is, co c ^ Acj>/At or A<j) = a) c At = 27rAt/T c ; so 
the definitions are consistent. The rest of this discussion will concern time 
uncertainty expressed as phase jitter. 

The first question regarding a telemetry system with coherent detection 
is: How much phase jitter can be tolerated before the system performance is 
degraded? Before one can answer that question the character of the phase 
jitter must be known. 


Characterization of Phase Jitter 


It should be emphasized at this point that discussing only M phase n 
jitter restricts the type of oscillators considered. In particular, the 
assumption is that the oscillator average frequency, denoted by 

^ - * c m * 

does not drift in the interval of interest (up to 100 s) . Thus the time 
waveform of the oscillator can be represented as 

G cos [2?rft + <j> (t) ] , | <\> (t) | <_ 7T rad (1) 


where G is an arbitrary constant, t is time, and <j>(t) is a realization of 
a random process (described as A<f> above) with an expected value E(<f>) = 0. 
The term <f>(t) accounts for the difference between the actual time of zero 
crossing of the waveform and the expected time of zero crossing of the func- 
tion cos 27rft. If the waveform _is phase modulated by a telemetry signal, it 
takes the form M(t) = G cos [2irft + <j>(t) + 0(t)] where 0 (t) is the telemetry 
information. A common way to characterize <j> is by its power spectral 
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density, which is useful (for our purpose) because in the phase domain the 
telemetry information can be separated from the noise if the power spectral 
densities of each are compared. Let S^m) be the power spectral density of 
ij>(t). Many researchers in the field of oscillator stability have studied 
S^fm). The resulting characterization has been reported in references 1 and 
2, which list numerous references, so the results will be stated here only 
briefly. The power spectral density of the phase jitter can be represented 
by the following polynomial: 


(m) = A/m 3 + B/m 2 C/a > + D 


( 2 ) 


in which A, B, C, and D are constants. The terms of the polynomial can be 
explained as follows. 

Within the tuned circuit of an oscillator, perturbations (noise) cause 
frequency changes. Two sources of noise are characterized. Because of ampli- 
fier gain variations and variations in the crystal itself there is a flicker 
(1/m) component of noise and, because of thermally induced voltages in the 
circuit losses, there is a flat component of noise. Together these components 
form t a frequency power spectral density S^fm) = A/m + B (rad 2 /s 2 )/rad. Since 
m = $, the equivalent phase power spectral density S<fi(m) can be derived from 
the frequency spectrum. The equivalence can be readily shown (although the 
method is not mathematically rigorous) using Fourier transform theory. The 
transform of 4>(t) is: 

3 [$00] = j«3[<Kt)] (3) 

but 2 13[4> Ct) ] | 2 can be termed the power spectral density, S^fm). Thus the power 
spectral density, (m) , of the internally generated noise is 

S* (<*0 =^S^(m) + ^rad 2 /rad (4) 


Outside the tuned circuit, amplifiers and multipliers operate on the 
oscillator signal, and perturbations in these circuits cause phase jitter. 
These perturbations can also be either flat or 1/m. This power spectral 
density is labeled S<j, o (m) = C/m + D rad 2 /rad. If <j>^ and <p Q are statisti- 
cally independent, the combined power spectral density is 

S<J>i + <p 0 (m) = A/m 3 + B/m 2 + C/m + D rad 2 /rad 

where the constants A, B, C, and D will be different for different circuits 
and also dependent on the choice of components. 

Equation (2) gives an idea of the qualitative character of the modulation 
phase jitter. For large frequencies, S^frn) appears to be flat and for 
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components near D.C. it appears to be steeply rising. The assignment of 
quantitative values of phase jitter to the oscillators and carrier in the 
design of a telemetry system can only be done by measurement. Unfortunately 
the implementation of such a measurement is not easily accomplished. 

The probability density function of phase jitter is gaussian with respect 
to some reference phase if the phase domain is unlimited. However, when 
dealing with periodic functions of phase, as in a phase tracking receiver, we 
think of phase as being restricted to the region from -it to +ir radians; 
thus the density must be modified by "folding over the tails" of the gaussian 
probability density function. Viterbi's formulation (ref. 3, p. 90) for the 
modified probability density function is 


a cos <}> 


= l + l * " <5) 

1 /• TT Qi COS 0 

where I is the zeroth-order modified Bessel function I Q (a) = — / e 

u TT " r\ 


de 


(see eq. 9.6.19, p. 376, ref. 4). In the ensuing discussion, a linear model 
will be assumed in which a = 1 /a| and is the variance of <j>. It should 

be realized that, for the operating range of interest for a phase tracking 
receiver (a^ < tt/4 is a reasonable bound for "in-lock" conditions) , this 
density looks much like a gaussian density (see figs. 2 and 3). 



Figure 2.- Comparison plot of two formulations of the probability density function of phase jitter; 
n(0) = ( e -^ 2 / 2 )/(27r/ a y / J , p(0) = (e« cos d>)/ 2 7rl 0 (a), a = l/a$. 
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Figure 3.- Comparison plot of two formulations of the probability density function of phase jitter; 
n(0) = (e^l^Kln/ayl 2 , p(0) =(e« cos *)/2irI 0 (a), a =l/aj. 


Total Phase Jitter in a Telemetry Receiver 


The discussion so far has concerned the character of phase jitter in an 
oscillator or transmitter; now it will be shown how this jitter is manifested 
when filtered by a phase lock loop receiver. First, however, the functioning 
of a phase lock loop (PLL) will be reviewed. Although the receiver is more 
complex, it will be sufficient to assume it to be a simple second-order loop. 
The PLL tracks, or follows, both the phase and the frequency of a signal and 
effectively filters out all but a narrow band of noise centered on the signal 
frequency. 

Three components are used to mechanize the PLL: A phase detector 

multiplies the input signal by a reference frequency to extract the relative 
phase error; a low pass filter removes all but the slowly varying or "D.C." 
components of the phase error; and a voltage controlled oscillator (VCO) 
generates the reference frequency for the phase detector. The VCO changes 
frequency in proportion to the output of the low-pass filter. A linear, phase 
domain model of the PLL will be used for analysis because it is directly appli- 
cable for small phase jitter, and the linear model is required for Fourier 
transform analysis. The phase-to-voltage transfer function of the phase 
detector is derived from the low-frequency part of the trigonometric identity 


sin(oi c t + cos (oo^t + <J> ) = 1/2 sin(<(K - <fr ) + 1/2 sin(2co c t + <Jk + <J> q ) 
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and from the approximation ((j)^ - fj> 0 ) small, sin(4>j - <j) Q ) - (<f>^ - <f> Q ) . The 
voltage-to-phase transfer function representing the VCO can be seen to include 
integration 1/jrn, since a voltage offset causes a frequency offset, which is 
a ramp in phase. Thus, we have the linear model: 


Adder Loop filler Units - 


v — i 

Integrator 

Sketch (a) 


M F (jw) h 


K d = V/rad 

F(jaj) = V/V 
K/joj = rad/s/V 


Define 


<f» 0 CjwD K 0 K D F(jaf) 

(j^O) jw + K Q K D FCjm) 


and for simplicity let = 1 . Then 



1 


H(jm) 


(6) 


(7) 


For a second-order loop, F(jco) = (jo)x 2 + l)/jwxi (which provides proportional 
plus integral control) , and the phase transfer function becomes 




H(jrn) 


j(OT 2 K/x i + K/x 1 
-u 2 + jmx^K/x^ + K/x ^ 

j(o26m n + w n 2 
-w 2 + jw26u) n + co n 2 


( 8 ) 


The phase error is 


1 - H(jco) 


to 2 - jo)26o) n - co n z 


( 9 ) 


where 6 = (x /2)(K/x )*/ 2 is the loop damping factor (6 = 0.7 for critical 

2. i / 2 

damping) and to n = (K/r^) 1 is the natural frequency of the loop. See 
figures 4 and 5. 
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Figure 4.- Phase transfer response for a high-gain, second-order loop where 5 = 0.707. 



Frequency, i*j/uj n 


Figure 5.- Phase error response for a high-gain, second-order loop where 5 = 0.707. 


To understand how these functions affect the total relative jitter in the 
receiver, one must differentiate between what we shall call internal and 
external sources of phase jitter on the incoming signal. "Internal' 1 means 
phase jitter that is truly part of the signal as described in the previous 
section, and for simplicity it includes any inherent VCO instability since we 
are concerned with relative jitter between signal and VCO. "External" jitter 
means phase jitter arising from thermal noise and is inversely proportional to 
the signal-to-noise ratio in the receiver PLL. The portion of these types of 
jitter which appear at the receiver phase detector will be designated <j>. 
and 4> x , respectively. 


The frequency spectrum S. is flat 

^X 

frequency is tracked out by the receiver, 
receiver phase detector due to this phase 


and the portion near the signal 

Thus the total power out of the 
jitter is proportional to 
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|H(ju>) | 2 du> 


( 10 ) 


or 



+* - ( 2 ’ b l ) s *, 5 7 (I), 


(ID 


where (N/S)^ is defined as the noise-to-signal ratio measured in the loop 
(ref. 5), is the standard definition of a single sided noise power band- 
width for a PLL, and is half the intermediate frequency bandwidth in 

rad/s. For 100 w n , and with a critically damped PLL, Bl = 0.53u) n Hz 

(ref. 5). (The bandwidth limitation is due to the intermediate frequency 

bandpass of the receiver which may be thought of as prefiltering the linear 
model of the PLL.) 

The frequency spectrum (oo) as described in the previous section is 

S (a)) = -^3 + + D 

(the C term is negligible in most oscillators) and the portion of this 
spectrum which appears as phase jitter in the receiver is that which is not 
tracked by the PLL (i.e., the phase error). Thus the total power out of the 
receiver phase detector due to this phase jitter is proportional to 


0)-j -p 

4>? n = f S + (<o) |l - H(jto) | 2 du> (12) 

o 

The exact computation of this integral is tedious, but the result can be 
bounded and the important terms can be found by approximating the phase error 
function from figure 4. Let | 1 - H(jtu) | - (o) 2 /co n z ) for a> <. a) n , and let 
|l - H(jo)) | - 1 for go > oo n . The approximation to the above integral is 
defined by U. z : 

J n T7 



U fn = / I 1 " H( »l 2 d(0 + V^l 1 - h O)I 2 d “ 

‘f 0 n ^ 4 D )(^7) 2,i “ 7 "(^ + ^ + D ) cl)2d “ 

U II w n 

i a) 100 %/ a R \ 

= 4 j* n (Au) + Bw 2 + Doj 4 )dco + J* \co^ + "^"2 + Djdw 

Wn *o *'o) n 

1 Taoj 2 Bw 3 Doj 5 1 w n r-A B .."l 10043 !! 

= ^L— 2 — *— 1 + b‘“ +D x 


r A B Du n 1 r A 

= 5 + + I + I ’9 - 

L 2 % 


i + r— 

J |_ 2a) n^ 


3u n 5 

Neglecting terms in A, B, and D of 1% or less, 


“n 


§ — 2 + ~ - - Doj„ + lOODm. 


2xlO H w n z w n 100(u n 


n] 


779 A A B B 

U-f-n = 0 2 + o 5 + 7 + + 100Do) n 

in 2(jo n z 2oa n ^ 3oi n oo n 11 


(13) 


(14) 


or 


2 = -At- + + lOODm, radians 2 


in oj n 2 3to n 


n 


Then 


U 


i 


n w C - m, ~ 2 + '7^ — + 100Doo n radians 
RMS w% 3u> n n 


(15) 


where 


U 


RMS " 



From measurements on the Pioneer spacecraft transmitter (see measurement 
results, p. 22) typical numbers can be assigned to the constants A, B, and D 
at S-band frequency to give a better perspective of the important terms for a 
narrow PLL (2Bp < 100 Hz). Typical values are A = 0.5 rad 5 /rad, B = 10 3 
rad 4 /rad, D = 2xl0~l° rad 2 /rad. Thus the A term clearly dominates for 
co n < 100 rad. Figure 6 shows that the approximation for |l - H(jto) | (A/w 3 ) is 
an upper bound for this term. It has been shown by approximation that only 


10 




Figure 6.- Power spectral density (PSD) for the relative phase error in a phase lock loop where the predominate 
phase jitter has PSD form A/co 3 . The PSD is in dB relative to A/c o n 3 . 



Figure 7.- Total RMS phase jitter in a phase lock loop where the PSD of 
the relative phase jitter is 0.5/co 3 rad 2 /rad. 


the A term is important, 
and fortunately its contri- 
bution to <f>RMS can be 
calculated accurately (see 
appendix A) . Figure 7 
shows Uj and <J> RM s for 
the A term versus 28^. 

The value of separating 
the two components of phase 
jitter lies in exposing 
conditions under which 
oscillator instability or 

cf)? n must be considered in 
the design of the threshold 
loop bandwidth of a 
receiver. Normally the 
receiver design threshold 
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is based on only the signal-to-noise ratio in the loop or <}>£. In most 
receivers the phase detector is preceded by a bandpass limiter which suppresses 
the signal and thereby reduces the loop bandwidth at signal-to-noise ratios 
near threshold. Therefore it is not surprising that <J>^ n can be a small phase 
jitter under strong signal conditions (where 2B^ is large) and still be large 
enough to reduce performance from the design value when the receiver is near 
threshold (where 213^ is small). Note the 10-times increase in <|>RMS f°r a 
10 to 1 change in 2 Bl in figure 7. 


Effect of Internal Phase Jitter on PCM Telemetry Performance 


In order to calculate the effect of internal phase jitter on telemetry 
performance, several assumptions are needed to further define the problem 
and limit the complexity of the solution: 


1. The receiver acts as a single conversion type with the PLL phase detector 
operating at the carrier frequency. The PLL is a second-order type 
critically damped as described in the previous section. 

2. The effect of receiver thermal noise in the tracking loop bandwidth will 
be ignored for this part of the discussion (i.e., only <t>^ n is 
considered) . 


3. No differentiation will be made between phase jitter originating in the 
transmitter and that originating in the receiver. All the jitter is 
shown in the receiver VCO. For the purpose of calculating the effect on 
telemetry error rate, only the relative jitter is of concern. 


4. The problem of subcarrier demodulation will be considered negligible; 
therefore the carrier input to the telemetry phase detector is modeled 
as being directly phase modulated with PCM. Also, perfect bit 
synchronization is assumed. The block diagram of the model is: 


G cos w, 




Loop I 
filter 1 


VCO 


sin b< + 't’inl t )] 

G cos[cu c t + M( t )] > (low poss filter [ -► e pl _ 


n( 1 ) - n, ( t ) cos io c t + n„( t ) sin <u c t 


Sketch (b) 


In the model, n(t) is the receiver additive 
noise represented by its in-phase and 
quadrature components, and epL is the 
baseband telemetry signal. Thus 


e Q = G cos [oo c t + M(t)] + n^(t) cos w^t + 
n q (t) sin u> c t (16) 


e = e sin [w t + <f>(t)l 

p O C 


(17) 


(For simplicity the subscript "in" is dropped from 4>(t).) Using trigonometric 
identities and neglecting higher frequency terms (the low pass filter removes 
all high order terms) : 


12 




+ 


sin <|>(t) 


cos <j>(t) (18) 


G n 4 (O 

e pL =ysin [<Kt) + M(t)] + — 2 — 


n a(t) 


but the first term can be written 


y cos 4>(t) sin M(t) + y sin 4>(t) cos M(t) 


Note : 

M(t) = ±d for PSK; cos ±d = cos d = y (y is a constant) 

Then 

r r n i (t) n q (t) 

e pL = y sin M (t) cos (j>(t) + y Y sin <j>(t) + — 2 — s ^ n ^(O + ~2 — cos 


In the above equation 

1 . y sin <j>(t) has a noise spectrum near DC, but in practice the data are 

on a subcarrier; thus this term is filtered out . 1 


2 . 


n (t) 

— sin <Ht) 

sketch (b) . 


n (t) 
+ _JL_ 


cos 4>(t) is 


an expansion of 


n(t) 

2 


as shown in 


Therefore 


e PL = T sin M (0 cos $(0 + ^y^ (19) 

where G sin M(t) = ±K; the magnitude of K is determined by the modulation 
index and the sign determined by the telemetry (+ corresponds to one; 

- corresponds to zero) . 

For correlation and matched filter detection of a binary signal such as 
this, the following decision variable is formed: 

q k =/ T [(±K cos 4>(t) + n ft) ) ] dt 

*o 

= ±K r T cos 4>(t)dt +f T n(t)dt (20) 

J o *0 

The symbol "one" is sent if is positive or "zero" if is negative. 

Now 

P e = Pr (error) where Pr( • ) is the probability of the event ( . ). 

x This is strictly true only for small <f>(t) where sin <J>(t) - <j>(t) . With 
a (j) = the sin <j>(t) will have energy at harmonics of the bandwidth of <|>(t), 

but only the higher harmonics would be in the bandpass of the subcarrier, and 
they are small enough to be neglected. 
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Then 


P = Pr(l) Pr(error/l) + Pr(0) Pr(error/0) 
e 

Also, by symmetry, 

Pr(error/l) = Pr(error/0) and Pr(l) + Pr(0) = 1 

thus 


P^ = Pr(error/l) 


Therefore it is necessary to use only 


T T 

q^ = K f cos <j>(t) dt + J* n(t)dt 

*o *o 


C21) 


To calculate the error probability, qi can be normalized by multiplying by 
1/KT without changing the decision point; then 


q = 



( 22 ) 


and 


P g = Pr(error/l) = J*° Pq(x)dx 


where p (x) is the probability density function of q. Unfortunately, 
finding the probability density of q is not easy. As can be seen from the 
above equation, q is the sum of two random variables which we will assume to 
be independent. 

1 f T 

The second term, -rpr; J n(t)dt, is the random variable derived from a 
1 l\ o 

white gaussian process; it also is gaussian with zero mean and variance 
o 2 =N 0 /2K 2 T. (N 0 = single-sided density; see appendix D for evaluation of this 
term. ) 


The probability density function of the sum can be found by convolution 
(see appendix E) : 


P 


q 



h (x-u)g (u)du 


(23) 


If 4>(t) is small (i.e., <)>(t) << 1) g(x) approaches a delta function at x = 1, 
then p (x) is a gaussian density with mean one. We know that for a gaussian 
Pq(x), P e is the error function for coherent PSK given by 
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I 



I 


p 

e 



1 


C\/5 TT 


e 


-Cx-n 

2o^~ 


2 

dx. 


N 


2K i 2 T 


By a change of variable, this can be written 


1 



(ref. 6) 


(24) 


Finding the probability density function of the first term is the main 
problem. The probability density function of ^ ^ cos <Kt)dt has not been 

found in general for the case where (w) is arbitrary and <j>(t) is assumed 

to be a gaussian process (see appendix B) . Therefore, the analysis will be 
restricted to the case when the symbol period T is short compared to the 
time required for <j>(t) to change appreciably. 


Symbol rates larger than the loop noise bandwidth are the main concern 
because the performance degradation for coded telemetry is of primary interest 
and the symbol rate for coding is relatively high. To illustrate this point, 
let P 0 = 0.001, corresponding to a normalized signal-to-noise ratio of 6.8 dB. 
If the loop bandwidth were equal to the symbol bandwidth and the power split 
equally between the carrier and the symbol, the loop signal-to-noise ratio 
would be at threshold (6 dB) which means, of course, that the receiver would 
hardly maintain lock. Actually lower symbol signal-to-noise ratios are used 
and less power is allowed for the carrier. Thus to maintain a higher signal- 
to-noise ratio in the loop, the loop bandwidth would always be a fraction of 
the symbol rate. 

Under these conditions cos (t) is practically constant over the 
integration interval T. Thus the following approximation holds (see 
appendix B) : 

i ,#T 1 T 

~ J cos <t> (t)dt - ■= cos <f>(t) / dt = cos <(> (t ) (25) 

1 o 1 *o 


and its probability density function is easily calculated (see appendix C) : 


where 


az 


g (cos <t>) = g(z) = 


TTl 0 (a) V7^ ' 


< 1 


(26) 


a 
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From physical reasoning g(z) is expected to be a function of the general 
shape shown in sketch (c) . The mean value of <f>(t) is zero, and the highest 
probability of occurs in the region about zero. Since cos is 

approximately 1 in this region, the highest probability of cos <f>(t) occurs 
near 1. The probability of realizing larger values of the magnitude of cf>(t) 
decreases smoothly according to the gaussian process, as the magnitude of <|>(t) 
increases. Therefore, the probability density of cos <{>(t) would be expected 
to decrease smoothly as |<f>(t)| increases. 



Sketch (c) 


This form for the phase jitter has been used to evaluate 


P e = y °h (x) *g (x) dx 


on a computer for various values of 


and % if = *RMS 


(See appendix E for the details.) Figure 8 illustrates the results. 



Three important conclusions can be derived from the results. First, if 
the total phase jitter is low C^RMS < 0-28 rad), the telemetry performance is 
nearly ideal over the range of P e of interest (P e > 5xl0~ 4 ). The performance 
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matches the theoretical error curve for coherent PSK within 0.2 dB for 
tf>RMs = 0.2 rad. Second, for higher values of total phase error, the normalized 
signal-to-noise ratio for the telemetry symbols rapidly loses influence on the 
error rates. It can be shown that for a given ^rms there is a P e which 
cannot be reduced by any increase in signal strength (see appendix E) . For 
example, if ^rms = 0*56, the minimum P e attainable would be approximately 
2.0><10" 2 . Third, in figure 6, which is drawn for oscillators such as those in 
Pioneer VI through IX, it should be observed that Aversely propor- 

tional to the receiver loop bandwidth in the range of interest (2B L <_ 12 Hz). 
Thus, for any given maximum allowable probability of error in telemetry, a 
minimum receiver loop bandwidth is established. For other oscillators and 
other ranges of P 0 and 26^, a proportional relationship between 2Bp and 
<t> ms will not be given in general; however, the same min-max relationship 
exists between 2Bp and P e . Fortunately, this bound often lies outside the 
range of P e and 2 Bl of interest to the designer. For Pioneer telemetry, 
the receivers currently used have a minimum loop bandwidth of 3 Hz and thus 
the effect of Pioneer oscillator instability on the telemetry is negligible 
as shown by the analysis results in figure 8 column B. 


THE OSCILLATOR INSTABILITY MEASUREMENT SYSTEM 


The difficulty in measuring <}>(t) is that one must generate a reference 
carrier equivalent to that of the test oscillator but without the <|>(t) phase 
jitter. As has been demonstrated above, in communications a time reference 
must be generated to measure or extract any phase modulation, including PM 
telemetry. In fact, the precision reference source built to extract <f>(t) 
contains a phase lock loop, just as a coherent PM receiver is basically a PLL. 
Of course, the design is different because the emphasis is on reducing the 
self-induced phase jitter to a minimum and on having a high sensitivity to 
phase variations. The resultant extracted phase jitter is the sum of the 
jitter from the two sources, and hopefully the phase jitter from the measuring 
instrument will be small. One such device was developed and built, under 
NASA-Ames contract, to measure phase noise in the telemetry link for the 
Pioneer series of spacecraft. It is called an Oscillator Instability 
Measurement System (OIMS), and consists of an extraction section 
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AM - 10%/V 



Figure 9.- OIMS extraction section. 


(see fig. 9) that extracts phase jitter from an intermediate carrier frequency 
of 10 MHz ±1 percent and a frequency conversion section that mixes various 
carriers, used in the Pioneer telemetry link, down to 10 MHz. To understand 
the design, one salient point must be demonstrated about phase jitter. If 
the passband for the carrier is always broad with respect to and the 

carrier frequency is multiplied by n, via a nonlinear device 

A cos n[2nft + <Ht)]= A cos [2imft+ncJ> (t)] nj) 


which means that the power spectral density is scaled up by n 2 but the 

envelope shape is unchanged. Likewise, if n is less than 1, say n = 1/D, 

is scaled down by 1/D 2 . For low phase jitter, the PLL for the OIMS has 
been built around a frequency standard with low phase jitter. Mixing its 
frequency with a relatively noisy synthesized frequency, which is partly 
voltage controlled and frequency divided by D = 100 (1/D 2 = 10~ 4 ) , produces 
a combined frequency which is variable. The divided synthesized component 
has a low phase jitter, and therefore the phase jitter is still small compared 
to most oscillators. 
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Figure 9 shows in detail how the extraction section of the OIMS functions. 
The 10 ±0.1 MHz reference is composed of two sources; a 9 MHz signal which has 
the stability of the frequency standard, and a 1 ±0.1 MHz signal which has 
components from the standard, the synthesizer, and a voltage controlled 
oscillator or VCO (which is physically part of the synthesizer). In this 
simple block diagram a linear model is used to describe the nonlinear elements 
that add and subtract, and scale frequencies up or down. Filtering and 
buffering requirements are not shown in the diagram. 

The 9 MHz signal is generated by simultaneously scaling the 5 MHz from 
the standard up to 10 MHz and down to 1 MHz, then taking the difference. The 
other component of the reference 1.0 ±0.1 MHz signal is generated by combining 
the 5 MHz standard, scaled up to 60 MHz, with the 40 ±10 MHz output of the 
frequency synthesizer. This produces a 100 MHz ±10 MHz signal which is scaled 
down to 1 ±0.1 MHz. The frequency of the synthesizer output can be selected 
manually over a range of ±10 MHz and has automatic tracking, via the PLL, of 
up to 10 KHz from the VCO in the kilohertz decade of the synthesizer. This 
range, of course, is scaled down to ±0.1 MHz and 100 Hz, respectively. 

The extraction section has three outputs which are limited to baseband 
frequencies. Each is buffered by low-pass filtering amplifiers (not shown) 
with variable gain , but only the most sensitive output transfer functions are 
indicated. The only output of concern to this discussion is labeled PM for the 
extracted "phase modulating" signal or equivalently phase jitter. The 
extractor also outputs the related "frequency modulating" or FM signal and a 
signal proportional to any incidental "amplitude modulation" (AM) . 

The 10 MHz reference, which is generated as described above, will track 
the phase and frequency of the signal under test without introducing much 
phase jitter of its own. Figure 10 shows the residual power spectral 
density of phase jitter for the OIMS extraction section when operated in a 
self-test or common mode condition. 2 Under this condition the frequency stan- 
dard also provides the signal to be tested and thus only the phase jitter due 
to the PLL amplifiers, phase detector, VCO, and synthesizer is displayed. 

Figure 11 shows the power spectral density of the phase jitter for a 
similar external frequency standard as extracted by the OIMS at 10 MHz. If 
one can neglect the residual, and both standards contribute equally to the 

2 Note that these plots, which were taken directly from low frequency 
spectrum analyzers, are in units of S^(f) rad 2 /Hz versus Hertz. The 
conversion to S^(io) is simple if one remembers that the power in a fixed 
bandwidth is the same, then 

Af c S < j,(f) = Au c s + (w) 

Af c S^(f) = A (2irf c ) S<j) (io) 


therefore 


S<j> (f ) = 2tt S^ (u3) 


(28) 


The crosshatched area indicates the band of uncertainty in the measurement. 
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Figure 10.- OIMS common mode power spectral density, S^f), of residual phase jitter. Spectral lines 
shown are harmonics of the AC power frequency. 


-80, 


I I i I I I I 1 | 1 1 — I Mll!| 1 1 1 I i I i I 


-90 


r— 


-lOOh 


2=. -120b 


-130b 


- 140 b 


-ISO' 


t L I I 1 I I I I l I l l I I l I I 


10 100 
Frequency, Hz 


J I I .1 .1 1.1 


Figure 1 1 OIMS extraction of S^f) for a frequency standard referred to 10 MHz. Spectral lines are 

omitted for clarity. 
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jitter, the amount associated with the OIMS standard is 3 dB less than this 
plot. Note that the phase jitter of the standard is greater than the common 
mode residual. 

The other section of the 0 IMS is the frequency conversion section (not 
shown) and its function is to heterodyne, up or down, signals under test into 
the range of the extractor (10 ±0.1 MHz). The phase jitter on the local 
oscillator in the converter also must be as low as possible since, when it is 
mixed with the signal under test, the resulting phase jitter power is the sum 
of that originating in the local oscillator and from the input signal. Thus 
all the local oscillator frequencies needed for tests are directly synthesized 
in the converter from the frequency standard used in the extractor section. 
Synthesis is performed in the same way as for the frequencies needed in the 
extractor, so the frequency converter is not described in detail. 


MEASUREMENT OF PIONEER OSCILLATION STABILITY 


Most measurements of the Pioneer oscillator stability were made at the 
intermediate frequency of 114.6 MHz. Thus a mixing or local oscillator (LO) 
frequency of 104.6 MHz was generated for this purpose in the frequency 
converter. The frequency multiplying and mixing from the standard's frequen- 
cies (5 MHz, 1 MHz, and 100 KHz) are represented by the following equation: 

LO = [(20x5) + (4x1) + (6x0.1)] MHz. The phase jitter on this local oscillator 
is dominated by the first term of 100 MHz. As was explained above, the phase 
jitter of this term has the same spectral shape as that of the frequency 
standard (fig. 11), but the magnitude is scaled up by n 2 = 10 2 or 20 dB. 

Figure 12 shows a plot of (f) vs f for the Pioneer oscillator taken 
directly from OIMS measurements of several oscillators of the same design; 
notice that the spectrum is considerably above S<j, for the local oscillator 
for frequencies below 100 Hz; thus the correct measurement of the Pioneer 
oscillator phase jitter is assured in this region. Also notice that the shape 
of S<j, for frequencies below 100 Hz is 1/f 3 . The rapid roll-off of for 

frequencies above 1.5 KHz is caused by narrow band-pass filtering of the out- 
put of the Pioneer oscillator. It is this plot that provides the basis for 
the terms A, B, and D used in the calculation of the total phase jitter in 
the Analysis section. 

From measurements made with the OIMS it was possible to simplify the 
general problem of calculating the effects of phase jitter on the telemetry 
performance by neglecting terms in the polynomial for the power spectral 
density. Knowing the actual magnitude and power spectral density of the tele- 
metry oscillator phase jitter made it possible not only to solve the specific 
problem of determining the oscillator's effect on Pioneer telemetry, but to 
draw general conclusions about conditions under which oscillator phase jitter 
is a critical parameter. 
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Figure 12.- OIMS conversion and extraction of S^f) for a Pioneer VI-IX type of transmitter oscillator 
multiplied in frequency to 1 14.6 MHz. Spectral lines are omitted for clarity. 


CONCLUDING REMARKS 


It has been shown quantitatively that short-term oscillator instability 
degrades telemetry performance. Instability was defined as a phase jitter 
with a general power spectral density, (go) = A/co 3 + B/oo 2 + C/co + D, and the 
oscillator instability measurement system for extracting the phase jitter of 
an oscillator was described. The instability in the Pioneer VI through IX 
telemetry transmitter circuits was used as an example, to derive an expression 
for the total phase jitter expected from a telemetry receiver as a function 
of receiver loop bandwidth. Finally it was predicted analytically that even 
under strong signal conditions, oscillator instability can degrade telemetry 
performance if the receiver loop bandwidth is narrow. Unfortunately the 
complete theory was not verified with data since in the Pioneer telemetry 
system the effect of phase jitter was negligible. 

The oscillator instability measurement system was key to these results. 
Measurement of the power spectral density of the phase jitter on the oscillator 
was a necessary first step. The problem of the effect of oscillator instabil- 
ity on telemetry is too broad a topic to treat generally, and data were needed 
to limit the parameters properly. The problem has been solved for a specific 
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case in great detail so that it will be easy to apply the results with 
confidence to another similar telemetry problem. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, March 11, 1971 
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APPENDIX A 


CALCULATION OF THE MEAN SQUARE VALUE OF THE 
INTERNAL PHASE JITTER 


In the Analysis section, the expression for the phase jitter in a phase 
lock receiver was derived. Presented here is the evaluation of the integral 
for the phase jitter, which is given by: 


<j»2 n = f lf (co) 1 1 - HCjuOpda) (Al) 

o 

The above expression represents the total phase jitter of <j>i n in the 
sum- of- squares sense. It also represents the vari ance of if $± n is a 

process with zero mean; a third interpretation of (j>? n would be that it is 
the modulation power loss on a carrier frequency due to oscillator phase 
jitter. The term S^(oj) is the single sided "power" spectral density of <f>i n 
and has units of radians 2 per radian. To relate S^fco) to a measurement made 
using the OIMS (see p. 18), it should be noted that S^fco) = (1/2tt) S^ff) where 
S^ff) has units of radians 2 per Hertz. Since (f ) is measured by a sweeping 

wave analyzer, which sweeps only positive frequencies, S^ff) is twice the 
density used when integrations are made over the whole frequency domain. 

To evaluate the integral, expand 1 - H(jco) and write the square of its 
magnitude 


1 - HfjoO 


! 2 _ 


or 


ju) 2 Sw n - to n 2 


(A2) 


Since 



Z 1 


Uil 

ll - 


z 2 


i z 2 i 1 

1 

and | z | 2 = zz* 

where 

z* is 

the 




|l - H(j 

») I 2 


12 - 


or 


L 2 ! 2 

joj 2 6(jd n - 0J n 


2 12 


(co z - co n z ) z + 46 2 co n 2 co 2 
.4 


oj 


(co 4 - 2o) n 2 o) 2 + co n 4 ) + 46 2 oo n " 


CO 
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But for a critically damped PLL, 26 = \!T, 46 2 = 2; thus 


|1 - HO «»| 2 ■ > (A3) 

The other factor in the integrand is which has been shown to have 

the general expression A/w 3 + B/to 2 + C/u> + D. However, approximations show 
that for loop bandwidths of interest, only the A-term in is 

significant; thus 






du 



(A4) 


Performing a change of variable, let z = co 2 ; then dz = 2 to dto and for 
0 <_ a) < u>j_f , 0 <_ Vz~ <_ or 0 < z < oj 2 ^ 
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2 J , 

c 

= A / 
? J , 


Jfi if 

2 03 do3 

•L 

Ll 4 
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“if 

dz 

Jq 

z 2 + 03^ 
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_A 
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P* tan 1 ~ Zi 2 ~\ 
L“n “n J 


03 ? ~ 

if 


2oj. 


A 03? r 
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— 2 tan 1 “ 


n -* o 


n 


03 . 


n 


(ref. 7, p. 294) (A5) 

(A 6 ) 


so that 


<P 


A 

RMS “ 



(A7) 


A 

Expression (A7) is exact for the total RMS relative phase error if ( 03 ) = 3 

in the region of interest. (However, the proper application of these results 
is predicated upon the PLL maintaining lock.) 
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APPENDIX B 


APPROXIMATION OF THE INTEGRAL OF A DECISION VARIABLE 


In the Analysis section, the decision variable q was used in the 
expression for the probability of error. This decision variable is the sura of 
two random components whose distribution must be determined. This appendix 
concerns the approximation 


T 


r T 

J cos $(t)dt - cos <j> 


(Bl) 


which will physically be justified for the choice of parameters, 
alternate approximation 

if f cos 4»Ct)dt - if fj [l - |(}) 2 Ct)]dt 

is not advantageous in this case. (The subscript "in" again has 
from $(t) for simplicity.) 


The 


(B2) 


been dropped 


The integral of a physically occurring random process, such as cos <j>(t), 
over definite limits of time, is a random variable (see appendix D) . Thus in 

1 f T 

the approximation J cos <j>(t)dt - cos <f>, the random variable 
1 rT 1 ° 

Y J Q cos cj)(t)dt is replaced by the random variable cos <j> that is obtained 
from the random process cos <J>(t) by fixing the time t = t Q (with 0 < t 0 <_ T) . 

1 rT 

Proceeding with the discussion of the quantity — cos (j>(t)dt, 

consider the following hypothetical experiment. Suppose that a particular 
time function (realization) from the random process cos <j>(t) is the input to 

an averaging device that performs the operation ^ / T . Thus the output is 
1 rT l o 

j J q cos <j>(t)dt, and the question arises. When is it nearly the same as the 

instantaneous value of the input? Clearly, the output will be almost the 
same as the input when the input does not change appreciably in the interval 
of time from 0 to T during which the device performs its averaging operation. 
Consequently, the relationship between the fluctuations that the input can 
undergo and the averaging time T needs to be examined. Since in reality T 
is the (telemetry) symbol period and it is fixed by the given bit rate it 

is preferable to study the changeability of the input. Furthermore, since 
the input cos 4>(t) to the hypothetical averaging device is a function of 
<f>(t) , it suffices to study the fluctuations that <f>(t) can undergo. Thus the 
power spectral density that contains information about the frequency 

content of <J>(t) will be studied. 


In the analysis section, it was explained that the phase jitter process 
<f>(t) in a phase lock loop receiver has the power spectral density 


s <t» 


A<a). 


n 


+ or 


(B3) 
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over the frequency range of interest, and S, is negligible outside of this 
range. Since decreases as 1/co 3 when w increases, <J>(t) does not con- 

tain any high frequencies within the significant range. This implies that 
cos <j>(t) varies only slowly so that its instantaneous value is close to the 
average value over a reasonably small interval of time. Defining the cutoff 
frequency f c (in Hertz) as the frequency at which the value of the power 
spectral density Sa is 3 dB below its maximum value, a quantitative relation- 
ship can be obtained between the length of time T over which the device 
averages the input and the frequency characteristics of <j>(t) . Since the aver- 
aging device gathers information about the input only during a T-second inter- 
val, an analogy to the Sampling Theorem for nonrandom signals may be drawn. 
According to this theorem, a time function containing no frequencies higher 
than f c Hertz can be completely characterized by values of the function at 
instants of time separated by l/2f c seconds. Applying the theorem it be- 
comes evident that the time function cos <Kt) does not change much during the 

1 T 

interval T seconds when T < l/2f c . Thus the output ^ f Q cos <J>(t)dt from 
the device performing the operation y jf^ will be almost the same as the 

input, cos <j>(t), for some fixed t between 0 and T < 1/2 f c . Since the time 
function cos <j>(t) does not fluctuate much during the T-second interval, it 
is clear that an arbitrary point t Q , with 0 < t Q <_ T, may be chosen to 

represent the average value ^ cos <J>(t)dt by the instantaneous value 

cos <|>(t 0 ) of the time function cos <f>(t). Writing cos <}> for cos <Kt 0 ) , the 
approximation 


1 

-y J cos <j>(t)dt - cos cj) 
o 

is obtained when T < l/2f c with f c the 3 dB cut-off frequency of Sy 
This approximation will be useful in appendix E where the distribution of the 
decision variable q will be determined. 

Now the disadvantages that arise from the alternate approximation will 
briefly be indicated (ref. 8, pp. 192-195) 




dt 


= 1 


1 _ 

2T 



First the approximation, cos cf>(t) - 1 - (l/2)<j) 2 (t) when <f>(t) is small, must be 

made; then one has to determine the distribution of the random variable 
IT 

2 j Jq 4> 2 (t)dt. The distribution of this random variable has been only approx- 
imately found when (i) (t) is a normal process and (i i ) for power spectra of 
the form S^(ai) = l/(a + cj 2 ) and of the form S<j,(cj) = constant for -cj q < cj < cj 0 , 
S<j>(u)) = 0 otherwise (refs. 9, 10). Since in this model the phase jitter 
process <p( t) is not of this type, rather 4> (t ) has Viterbi's distribution (as 
stated in the analysis section) and S^fu) = Acj/(cj n l+ + cj 4 ) , the wealth of 
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existing knowledge does not help. Fortunately the physical situation 
suggested a simpler approximation so that there was no need to use a number 
of approximations that do not apply directly to the model. 
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APPENDIX C 


DERIVATION OF THE PROBABILITY DENSITY FUNCTION OF cos <Kt) 


In this appendix the probability density function g of the random 
variable Z(t) ^ cos 4>(t), at a fixed time t, is derived from the probability 
density function p of the random variable <f>(t). Since <f>(t) represents the 
phase jitter process, it was reasonable to assume in the analysis that its 
probability density function has Viterbi's form (ref. 3) 


/ a cos <t> 
e 


p(<f>) =' 


2TrI 0 (a) 


-7T < (j) < TT 


otherwise 


(Cl) 


The real number <f> stands for one value of the random process <}>(t) at time t. 
Let z stand for the corresponding value of the transformed random process 
Z(t) = cos <j>(t) at the same time t. 

Now the probability density function g(z) will be expressed in terms of 
the probability density function p(<f>). Since z is the argument of the 
desired function g(z) and since the argument of the given function p((j>) is 
<j>, first <(> is obtained as a function of z so that p(<|>) can be written in 
the form p[<J>(z)]. From z = cos <j> it follows that <j> = <J)(z) = cos“ 1 z. Now 
for any given real number z Q that lies between -1 and 1, the equation 
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Figure 13.- The function cos _1 z. 


(C2) 


has two solutions, denoted by 
<J> = <J> 1 and <J> ~ (j) 2j which lie 
between -tt and tt, as shown in 
figure 13 for one particular 
choice of z 0 . This figure 
shows that the function cos _1 z 
is symmetrical with respect to 
the z-axis. This means that 
one of the solutions of the 
equation (p = cos" 1 z 0 is the 
negative of the other solution; 
that is. 


<f>2 = -<h 

which will subsequently be 
used in the argument of the 
function p (<f>) . 
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The transformation cos <}> = z of the value <f> (of the random variable 
c|)(t) for fixed t) to the corresponding value z (of the new random variable 
Z) transforms the probability element [p(<£i) + p (<j>2) ]d<j> to the corresponding 
probability element g(z Q )dz. In fact, the two probabilities are equal; 
that is. 


g(z 0 )dz = [p(<h) + p(<t>2)] d< t> 


(C3) 


because probability must be conserved under the transformation. This trans- 
formation of probability can be represented graphically as a transformation of 
area elements in figure 14. (Note that the sum of the areas of the shaded 
elements equals the area of the crosshatched element.) 



Figure 14.- The transformation of probabilities. 


As can be seen in figure 14, the magnitude of the derivative 
d<j>/dz = (d cos -1 z)/dz at the point z G i s sc ale factor by which the 

shaded areas are transformed to the crosshatched area. From this geometrical 
reasoning, which is based on the principle of conservation of probability, the 
general relationship follows 


g(z) 


dj> 

dz 


[pOi) 


+ PO2)] 


(C4) 
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which expresses the transformation of probabilities under the transformation 
cos 4> = Z of random variables. Computing 


dd) 


Id cos^z 

dz 


dz 


1 

Vl - z 2 


z ^ ±1 


(C5) 


and realizing that p is an even function, that is. 


P0l) = P(-<I>2) 
= P(<t>2) 


(because <j>2 = the formula 


g(z) 


d<t> 

dz 


2p (<)>) 


1 

Trl 0 (a) 



(C6) 


is obtained. The variable $ is restricted to the interval -tt < <(> < tt, 
and z = cos (f> is restricted to the interval -1 < z < 1. When either 
z<~l or z > 1, there are no values of <j) corresponding to these values 
of z. Thus the probability that either the random variable Z < -1 or Z > 1 
is 0. 


This completes 
probability density 


the proof that the random variable 
function 


Z = cos 


has the 


g(z) 


1 

Trl 0 (a) 



V 0 


-1 < z < 1 


otherwise 
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Figures 15 and 16 show this probability density function for various values 
of the parameter a <j, = l/\Za] and figure 15 also shows the transformation of 
PO) to g(z) . 



Figure 15.- Transformation of PDF under the change of 
variables from “0” to “cos 0.“ The PDF’s are drawn 
to scale for = 0.56, a = 0.312 (<y = 1/a^); 
p(0) is on the cylinder of unit radius, perpendicular to 
the x-y plane; 0 is on the unit circle; g(cos 0) is in the 
x-z plane; cos 0 is on the x-axis; dF denotes elements 
of probability. 



Figure 16.- Probability density function; 
r(z) = e a7 /7Tl(a)( 1 - /. 2 )' n , a = 1 \a \ , 
for two values of o^. 



APPENDIX D 


DERIVATION OF THE PROBABILITY DENSITY FUNCTION OF J_ f T n(t)dt 

KT *o 

The decision variable q was used in the expression for the probability 
of error associated with a telemetry signal. This decision variable is the 

1 /*T 

sum of two random components. The first random component ■=? J cos <j>(t)dt 

1 O rp 

was treated in appendixes B and C. The second random component f n(t)dt, 
where K is a real constant, will be discussed here. 0 

It will be necessary first to give a physical meaning to the quantity 

1 T 

rr r 1 n(t)dt. The random process n(t) consists of a collection of determinis- 

tic time functions (realizations) that can be represented in the set 
{ni(t), n 2 (t), . . . , n^(t), . . . , nk(t), . . . }. In any one particular 
experiment, one time function, n^(t) is integrated over an interval of T 

seconds yielding the number / T nj[(t)dt. Then this number is divided by the 

O 1 j»T 

integration time T giving the number ~ / nq(t)dt which is recognized as 

* o 

the time average of the function n^(t) over the interval of T seconds. 

When this experiment is repeated, a different time function, n^ft) is involved 

1 m T 

whose time average ^ J Q %(t)dt computed. Continuing the experiments in 

this manner, a collection of time averages is obtained which may be repre- 

( 1 /* T 1 /*T 1 T 

— J ni(t)dt, ^ j ri 2 (t)dt, . . . , i f l ni(t)dt, . . . , 

I j O 1 O AO 

Y J Q 1 n] c (t)dt, . . .1 . The elements of this set are real numbers that 

resulted from having performed the experiments 1,2,. ..,i,..., 

k Before each of these experiments it was not known which one of the 

time function averages would be computed. Consequently, it was also not known 
which one of the elements in the above set would be the outcome of the partic- 
ular experiment. This ignorance renders the set above a random collection, 

and it is represented by the symbol n(t)dt that is called a random 

variable in probability theory. Thus, for any random process such as n(t), 

from which the integral over a constant interval T of time is computed, such 
i* t 

as n(t)dt, a random variable is obtained. Figure 17 summarizes graphically 

these introductory remarks. The sum of the shaded areas above the t-axis minus 
the sum of the crosshatched areas below the t-axis represents T times the 
time average of the respective time function. Figure 17 is drawn for an 
arbitrary random process n(t); in the following discussion n(t) will stand 
for a white, Gaussian noise process. 

Now some qualitative statements follow about the random variable 


Y =T?F / T ”Ct)dt 
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Figure 17.- Typical time functions from the random process n(t). 


In the analysis it was 
assumed that the noise 
process n(t) obeys the 
normal probability law 
with zero mean and vari- 
ance a 2 so that the 
probability density 
function of n(t) is 



a V2 it~ 


for < n < This 
means that the noise 
process n(t) takes a 
value n k in the interval 
from a to a + dn with 
probability 


h(n k )dn 


1 


a \/5 t 7 



dn 


When the value n k approaches infinity, the probability lirn^ h(n k )dn of 

~n 2 *^k 

this happening tends to 0 at the rate e Thus it follows that the 

integral j*^n(t)dt exists and converges for almost all time functions in the 

set {ni (t) , n 2 (t), . . . } = n(t). Davenport and Root (ref. 11) show that 
under these conditions the integral of a normal random process is itself 
normally distributed. Since a normal random variable is statistically com- 
pletely specified by giving its mean and variance, it is necessary to deter- 
mine only those parameters in order to write its probability density function. 

2 The following is the determination of the mean E{Y}, and the variance 
o y of the random variable 


Y = KT f Q n ^ dt 

that represents 1/K times the time average of the random process n(t) over 
the interval of length T seconds. Since the expectation E is linear 
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E{Y } = E 



n(t)dt 



E{n(t) }dt 


Since the mean E{n(t)} of the noise process n(t) was assumed to be 0, it 
follows that 


E{Y } = 0 

In such a case the variance ay = E{ [Y - E{Y}] 2 } reduces to the mean square 
value E{Y 2 } which will be computed next: 

e{y2> ' E {(iSf / 0 Tntt)dt ) 2 } 

n(t)n(u)dt du} 


Since the double integral in the second equation above can be regarded as the 
limit of a double sum, and since a double sum is simply one large sum, the 
linearity of the expectation E allows us to interchange the order of 
expectation and integration. Thus 

E { (£ n(t)n(u)dt duj = E{n(t)n(u) }dt du 

Y y 

= f f R n (t,u)dt du 
O O 

where R n is the autocorrelation function of the noise process n(t) which is 
also assumed to be stationary so that R n (t,u) = R n (t-u). Since 
S n (f) = (N 0 /2) W/Hz for -» < f < “, it follows from the Wiener-Khintchine 
Theorem 1 that 



R n(t - u) = 3 _1 {S n (f) } 

N n 

= 6 (t - u) watts 


^he Wiener-Khintchine Theorem states that the autocorrelation function 
and power spectral density of a wide sense stationary process are a Fourier 
transform pair. 
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where the symbol 3 1 denotes the inverse Fourier transformation. Using one 
of the properties of the Dirac delta function, that is, 

f 6(w)dw - 1 


and with the limits on the integral suitably chosen, the following equations 
can be written: 



-T r T N q 


- u)dt du = J J 6(t - u)dt du 


Nq f 1 
2 J ~ 


du 


NqT 

2 


Thus the mean square value 


o N 0 

EtY * " 2K 2 T 


o 

is obtained which equals the variance ay of the random variable 



The probability density function of Y then is 


h(y) 


V (N 0 /2K 2 T 


-y 

2 (N d /2K 2 T) 


f or -oo < y < oo . This completes the determination of the probability density 

1 fT 

function of the random variable ~ J n(t)dt, taking real values y, that is 

o 

the second random component in the decision variable q. 
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APPENDIX E 


EVALUATION OF SYMBOL ERROR PROBABILITY 


In the analysis the expression 


P e = f ° PqOOdx (El) 

represented the probability of error where pq is the probability density 
function of the decision variable q. Since the random variable q is the 
sum of two independent random variables with probability density functions h 
and g, respectively, it follows from the Convolution Theorem 1 that 


PqO) = h(x)*g(x) 



h(x - u)g(u)du 


(E2) 


In appendix D the probability density function 


h(x) 


1 e 2(N D /2K 2 T) 

V(N 0 /2K 2 T)^27 


-oo < x < 00 (E3) 


and in appendix C the probability density function 

ax 

( g — -- , -1 < x < 1 

I ttI 0 (a) Vl - x 2 

g(x) =i (E4) 

I 0 , otherwise 

were derived. Since the function g(u) in the integrand of (E2) is zero for 
values of the argument u that are outside of the interval from -1 to 1 and 
since the argument x - u of h is unrestricted, the range of integration has 
only to extend from q = -1 to U 2 = 1. (This observation can be proved 
rigorously by the method outlined in the article cited in the footnote.) 


^he concept of convolution and its physical meaning are lucidly 
presented in the article by T. J. Healy, reference 12. 


38 



Ill 


Substituting the integral for Pq into the expression for the probability 
of error (El) results in 


P 


e 



h(x - u)g(u)du 


i 


dx 


(E5) 


Now P e will be expressed in terms of tabulated functions. Since the 
integrand h(x - u)g(u) is (absolutely) continuous for -1 < u < 1 and 
-°° < x < 0, the order of integration may be interchanged so that 


^e 




g (u) du 


(E6) 


Now define 


\/2(N 0 /2K 2 T) 


where u is fixed, so that 


dx = - V2(N 0 /2K 2 T)dv 


The limits of integration x = -°° and x = 0 become v = +°° and 


u A 


V 2(N 0 /2K 2 T) V £ 


Then the integral in brackets of (E6) is written 


/: 


h (x - u)dx 


i r°° 

f— •'vo 
VTT 


“ V A 

e dv 


(E7) 


having used the formula 


f f(y)dy =J f(y)dy 


The co-error function, erfc, is defined by (ref. 4, entry 7.1.2) 


erfc(v £ ) = 


oo 2 

-V^ , 

e dv 
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III 


so that 


1 /•» _v 2 1 

-J=£ e dv- — erfcfv ) 

The final expression for the probability of error is 


where 


and 


2 J- 1 [ V2(N 0 /2K 2 T)] 


g (u)du 


(E8) 


,au 


g(u) = 


irl n (a) Vl - u 2 ’ 


-1 < u < 1 


otherwise 


, . A 1 Z* 77 a cos <)> 

'• W = 7 Jo e 


dcj) 


is the modified Bessel function of order zero (ref. 4, entry 9.6.19). Also, 


a 


E{ *in> 


since the mean of f is zero. 

T m 



The final question concerns the form that the expression for P e assumes 
under strong signal conditions. Under such conditions the random variable 
Y = ^ jf T n(t)dt is much more likely to take on small values than large. 

Since a ^ = N 0 /2K 2 T the variance of Y (or the noise-to-signal ratio), a 

strong signal condition is then equivalent to the requirement that ay 
become vanishingly small. In such case the probability density h(y) of Y 
will be large about its mean E{Y} = 0 and will decrease rapidly for values 
y away from the mean. In the limit as ay approaches 0, all probability 
is concentrated at the mean, that is. 


2 

lim h(y) = lim - — e 2 °Y 

<7^*0 CTy-^O °y\ / 5tt" 

= fi(y) 
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From the sifting property 



u)g(u)du = g(x) 


it follows that 


PqOO = h(x)*g(x) = g(x) 


when a -*■ 0. Thus 


Pe 




g (x) dx 


because g(x) = 0 for x < -1. Substitution of the probability density 
function g yields the expression 


Pe 


1 /• o 

ffI oC°0jLi 


e ax 

Vl - x 2 dx 


for the probability of error under strong signal conditions. In this case P e 
equals the probability of cos </> taking negative values which is a function 
of the variance. Evaluation by computer showed that the previous expression 
for P e (containing the co-error function erfc) yields the same values for 
P e when ay is small as does the preceding expression. Figure 8 shows P e 
as a function of the standard deviation of (pj n and as a function of 

I £ ST (expressed in dB) . 

2(N q /2K 2 T) N q 
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